Singular elastic strains and magnetoconductance of suspended graphene 
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Graphene membranes suspended off electric contacts or other rigid supports are prone to elastic 
strain, which is concentrated at the edges and corners of the samples. Such a strain leads to an 
algebraically varying effective magnetic field that can reach a few Tesla in sub-micron wide flakes. 
In the quantum Hall regime the interplay of the effective and the physical magnetic fields causes 
backscattering of the chiral edge channels, which can destroy the quantized conductance plateaus. 

PACS numbers: 73.61.Wp, 73.43.Cd, 73.23.Ad 



The isolation of graphene monolayers [U [5] and the 
observation of the integer quantum Hall effect (QHE) in 
these systems [31 0] have made graphene a very active 
research topic O [6]. The integer QHE in graphene is 
remarkably robust and can be seen even at room tem- 
perature [T . The traditional setup to measure the quan- 
tized Hall plateaus involves at least four contacts: source 
and drain for the current, and two (or more) side con- 
tacts for the voltage. This scheme helps to eliminate the 
spurious contact resistance. The much higher mobility 
of suspended samples [5HTU] raised hopes for observing 
also the fractional QHE in graphene. Yet demonstrating 
even the integer QHE in a four-contact setup proved to be 
difficult in such samples. Only recently the integer QHE 
has been confirmed in suspended graphene [9j [HI [12] by 
reverting to a two-contact scheme. (Similar observations 
have been made for bilayers [13].) The fractional QHE 
has also been reported in these experiments [TTl [T^] . An 
artifact of the two-contact setup is the suppression of the 
quantized conductance [II1[T2], familiar from the QHE in 
semiconductors |14H17j . The reason why the nominally 
superior multi-contact scheme is less successful in sus- 
pended graphene has not been fully clarified, except for 
one theory that in small samples the side contacts had 
to be placed too close to the source and drain, causing 
admixture of the longitudinal and Hall conductances . 

In this paper we consider a different effect, which may 
also contribute to the lack of quantization: when the 
graphene sheet is under tension, the side contacts in- 
duce a long-range elastic deformation which acts as a 
pseudomagnetic field M{x,y) for its massless charge car- 
riers [5 , 18, 19,. The tension can be generated either by 
the electrostatic force of the underlying gate [201 [21] , by 
interaction of graphene with the side walls [52], or as a 
result of thermal expansion [531 IMl- Our main results 
are as follows. We show that B is concentrated near the 
corners of the contacts where it exhibits power-law singu- 
larities. It decays into the interior of the sample but for a 
reasonable 0.1% average strain in a 200-nm wide strip, B 
can remain of the order of a Tesla across its entire width. 
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FIG. 1. Conductance in units of Go = 2e^/h as a function 
of the filling factor v for samples of electron density 2.8 x 
10" cm"^ and widths W = 100, 200 and 400 nm. The side 
contacts are assumed to have width 21 = 30 nm and to pull 
laterally outward with the force P — 78 nN (co — 1) each. 

This leads to backscattering of the QHE edge states when 
the real magnetic field S is in a similar range, causing 
the erosion of the quantized conductance plateaus. We 
give an analytical argument that predicts that the QHE 
plateaus are destroyed above a threshold Landau level 
index Nc and the corresponding filling factor i^c'- 

N.^Omco'^^ikpWf/^ , z/, = 4iV, + 2, (1) 

where kp is the Fermi momentum in zero magnetic field, 
W is the width of the sample, and cq ^ 1 (see below). We 
also compute the conductance numerically. The results 
for kp = 0.94 X lO^cm"^ and three representative W's 
are shown in Fig. [I] For these three traces Eq. ([T]) gives 
(top to bottom) i>c ~ 38, 16, and 7.6, in agreement with 
the simulations. Based on these results, we suggest that 
observation of the QHE in suspended graphene with side 
contacts requires samples of width W > 200 nm. On the 
other hand, in small samples one can envision purposely 
using the pseudomagnetic field as a new tool for tuning 
electronic properties of graphene nanostructures [25-28] . 

We model a suspended graphene sample as an elastic 
membrane occupying the rectangle |x| < L/2, \y\ < W/2. 



The membrane has the two-dimensional (2D) Young 
modulus C = 340 N/m « 2100eV/nm^, the Poisson ra- 
tio a — 0.15, and the shear modulus — C/[2(l -f a)] w 
900cV/nm^. This membrane is supported by four con- 
tacts that have the same uniform height. We assume 
that the sample is in the state of the plain stress, i.e., 
we ignore the possibility of spontaneous wrinkling |29j . 
We model the two side contacts as rigid stamps of width 
21 < L centered at the symmetry axis a; = 0, see Fig. [2] 
The remaining parts of the y = ±W/2 boundaries are 
assumed to be free of tractions. Each of the stamps pulls 
graphene normally outward with the total force P. Ini- 
tially we will assume the stamps are rigidly clamped to 
the sheet. Later, we will consider other boundary con- 
ditions. At the X = ±L/2 sides (the source and drain 
contacts) the sheet is clamped: u ~ v — 0. Here u{x,y) 
and v{x,y) are the elastic deformations in the x and y- 
directions, respectively. At P = the membrane is sup- 
posed to be flat and unstressed. We are interested in the 
deformation that develops if the force P is finite. 

Our program is as follows. First, we obtain an ap- 
proximate solution of the posed elasticity theory prob- 
lem. Next, we determine the distribution of the effective 
vector potential according to the formulas [El [30l [31] 

K,^ I3'^{d.,u-dyv), Ky^-l3^dyu + d^v), (2) 

where = hc/e is the flux quantum (reduced by 27r), 
a — 0.14 nm is the separation of the nearest carbon atoms 
and P — d\n.t/d\na is the logarithmic derivative of the 
corresponding hopping integral t « 3 eV. 

ThetwoEqs. ^ can be combined using the complex- 
variable notation: 

A = A^. + iky = {2l3ip/a)dJJ , U = u + iv , (3) 
M = dxAy-dyA^^2lmd,_A, z = x + iy , (4) 

where the bar denotes complex conjugation. The last 
equation gives the pseudomagnetic field. As a final step, 
we calculate numerically the ballistic conductance be- 
tween the source and drain contacts assuming the system 
is subject to the total magnetic field Btot = B + M. 

Our elasticity theory problem does not seem to have an 
analytic solution but we can construct an approximate 
one from those for a half-plane. Consider therefore a 
membrane that occupies the entire lower half-plane y < 
0. Let its boundary be unloaded except for the interval 
—l<x<lm contact with a rigid stamp. The following 
general representation [3? is valid: 



2^U{z) ^ x(l){z)+(j){z) + {z-z)(t)'{z) , 
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where function 0(2:) and its derivative <&(z) = '\>' [z] are 
regular in the complex plane of z with a branch cut 
{—1,1). These functions are determined by the shape of 
the stamp and the forces it exerts on the membrane. For 
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FIG. 2. Schematics of the system. The rectangles at the ex- 
tremes of the X and y ax.es are the contacts, which are attached 
to a graphene sheet (the central object). The deformation of 
the sheet is strongly exaggerated. The variable shading and 
honeycomb lattice are added for aesthetic purposes. 



a stamp whose base is straight, parallel to the x-axis, and 
bonded to the membrane we have [32, 
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with TO = (ln>f)/(27r) w 0.14. The asymptotic behavior 
of '^{z) at large z — re*^ is given by 



$ z ~ -2zco^ -e 
p r 



r:>l. 



(7) 



same as in the Flamant problem: a point force applied 
to a 2D elastic sheet [33]. Here we introduced the di- 
mensionless coefficient cq = /3P/(47r/za) , for convenience. 
Substituting Eq. ^ into Eq. ^ and then to Eq. Q, we 
find that B decays as the square of the distance and has 
a peculiar angular dependence: 



S cos 46* -4 cos 26*) , r>/. (8) 



On the other hand, near the corners of the stamp, z — ±1, 
the solution Q is characterized by divergent oscillations. 
As discussed in Ref. [32], such oscillations are unphysi- 
cal. The linear elasticity theory usually fails at distances 
where they supposedly occur. Hence we consider the fol- 
lowing alternative to Eq. ([6]): 
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The corresponding elastic deformation is given by 



a / - HZ 2z — z 



const. (10) 



This happens to be the exact solution for a frictionless 
stamp with a rounded base [32] . It is physically relevant, 
free of divergences, and has the same universal far-field 
behavior, Eq ([t]). Since it is the far-field behavior that 
is important for the backscattering of the edge states in 



wide samples, we adopt Eq. ( 10 1 as our basic building 
block. (The subscript "1" in Ui is to remind us that it 
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FIG. 3. (Color online), (a) Pseudomagnetic field, in Tesla, 
induced by two stamps with parameters 21 = 30 nm, cq = 1 in 
a sample of width 200 nm. (b) Effect of this pseudomagnetic 
field on the edge states of the N = 2 Landau level at B = IT. 
The edge states are depicted as light gray ribbons of thickness 
£ ~ 26 nm. The arrows indicate their propagation direction. 



is for a single stamp.) We now construct the solution for 
the original problem simply as the sum 



U{z) = Ui{z~ iW/2) - Ui{~z - iW/2) . 



(11) 



Let us now estimate the coefficient cq in Eq. (10 1. 
Up to logarithmic factors, cq ^ /SeyyW/a, where Cyy — 
\U{iW) — U{0)\/W is the average strain. We see that the 
estimate [20] of Cyy ~ lO""* for a sample of size W ^ 1 ^m 
yields Cq ~ 1. In this paper we are interested in samples 
of smaller size, W ^ 100 nm. However, we think that 
Co = 1 is still a reasonable number if the pulling force P 
is determined by the direct interaction of graphene with 
the contacts, not via the intermediary of the gate un- 
derneath. The distribution of the pseudomagnetic field 
computed for such Cq and W using Eqs. (If, ^, ([To]), and 
(11) is shown in Fig. |3] As one can see, B ~ IT. This 
suggests that strain can strongly affect magnetotransport 
if B is comparable or smaller than IT. A more precise 
criterion can be derived as follows. 

The suppression of the A'^th quantized plateau (cen- 
tered at filling factor 4N + 2) is due to backscatter- 
ing between the counter-propagating edge states of the 
A^th Landau level. For A^ ^ 1 a semiclassical descrip- 
tion is valid, in which the edge states are visualized as 
channels of width ~ ^ = ^pjB that follow the con- 
tours £Ar = hvF\/2N \Btot{x, y)\ /f — Ep, where vp 
is the Fermi velocity and Ep — hvpkp is the Fermi 
energy. In the middle of the A^th quantized plateau 
Ep = hvpyJ{2N + l)B/(p , so that the edge states fol- 



lix,y) = B/{2N) = const. 



(12) 



These contours eventually terminate at the sample 
boundaries because B decreases away from the stamps. 
Upon reaching these termination points, the edge states 
continue along the boundaries, as usual [Fig. [sjjb)]. 

On the other hand, near the stamps the edge states 
veer into the bulk. The counter-propagating edge states 
can approach each other near the saddle points oiM(x,y). 
Using Eq. ^ one can show that for 21 <^W the saddle- 
points closest to the origin are located at y = and 
X = ±xo, where xq = W{1 - 2/V5)^/^/2 « 0.16MK. Near 
X = xq we find 

M{x,y)~-cocMx-xo)y/W\ ci « 2988 . (13) 

As the real magnetic field decreases, strong backscatter- 
ing and therefore the complete destruction of the QHE is 
possible if the distance between the counter-propagating 
edge states at the saddle-point becomes ~ £. Using 
Eqs. ([12]) andJlS] 
/cqCiN 
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this condition can be 
Substituting £ = y/2N ■ 



written as 
we 



1/kp 

recover Eq. (|I]). 

Now we discuss our numerical results for the ballis- 
tic conductance G shown in Fig. [l] The calculation was 
done by means of the recursive Green's function tech- 
nique [3S]. The system was treated as a tight-binding 
model on a honeycomb lattice with nearest-neighbor hop- 
ping t = ioexp(-/3|(5a|/a). Here [3 = 3.37 (US]) and 
Sa{z) = U{z + t) — U{z) is the complex-valued distor- 
tion of the bond that connects sites z and z -I- r in the 
unstrained lattice. Although we have referred to the side 
stamps as the contacts, in our calculations they are elec- 
trically isolated from the system, i.e., G is the conduc- 
tance between the source and the drain. Still, our results 
are representative of the four-contact conductance be- 
cause we made sure that contact resistances of the source 
and the drain are small. In the L oo limit where the 
source and drain are infinitely far from the strained re- 
gion these resistances vanish. We approached this limit 
by increasing L while keeping I and W constant until 
no significant L-dependence of G was seen [35]. As one 
can see from Fig. [T] the calculations show the expected 
erosion of the QHE plateaus as v increases and confirm 
the analytical estimate of the threshold filling factor i/c, 
Eq. ^. 

In conclusion, we want to mention a few caveats: 

(i) besides the pseudo vector potential A, the strain 
also creates a scalar one [TS]. However due to screening, 
it is likely to be less important [201 HZ] • 

(ii) We have assumed that electron transport is co- 
herent. If dissipation is present, then the two-terminal 
conductance is better understood in terms of mixing of 
the Hall and longitudinal conductivities [TTJ [Tni [35] ■ This 
effect is strongly enhanced in long samples |16j . 
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(iii) The point where a contact ends and a free edge 
begins can be considered a corner of angle a = tt at which 
the boundary conditions suddenly change from clamped 
to free ones. Divergence of strain near such corners is well 
known. The exponent A of the corresponding power law 
depends on the Poisson ratio a and the corner angle a. 
For a = TT and a = 0.15, Eq. ^ implies A = l/2± 0.14i. 
If we had instead a = tt/2, it can be shown [321 ES] that 
A would be 0.16. In this case, the strain at a distance of 
order £ from the corner is enhanced with respect to its 
value near the center of the flake by a factor of {W/£)^. 
For i ^ 10 nm and W ^ 1 fiui this factor is ~ 1.5, which is 
a sizable effect. Accordingly, the recently proposed [25- 
[28] "strain engineering" of graphene electronic devices 
must either take into account or take advantage of the 
strain singularities we discuss here. 
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